Molecular transitions in Fermi condensates 
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We discuss the transition of fermion systems to a condensate of Bose dimers, when the interaction 
is varied by use of a Feshbach resonance. We argue that there is an intermediate phase between the 
superfluid Fermi gas and the Bose condensate of molecules, consisting of extended dimers. 
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The remarkable experimental advances in atomic trap 
physics now permits one to study Fermi gases of atoms 
over a broad range of couplings from weak to strong. 
In particular, two-component Fermi systems have been 
made and studied using the two alkali metal isotopes 
40 K 0, H Q and 6 Li [J H @. In both cases, the Fes- 
hbach resonance phenomenon has been exploited to vary 
the coupling strength of the interaction between atoms 
and study the consequent many-particle properties. In 
both cases there is a closed channel near threshold who 
energy can be tuned with respect to the open two-body 
channels. A particularly interesting feature is the forma- 
tion of molecules or dimers as the effective interaction 
strength is increased. There are many questions than 
can be posed about a coexisting phase of bound pairs 
and that we would like a theory to address. 

Here we want to argue that the theoretical situation 
is rendered more complex by the fact that there are fact 
two kinds of pair states that must be treated in any the- 
ory that applies to the full range of couplings. These two 
kinds of pairs are spatially different and have different 
quantum numbers as well. The first bound state appears 
when the inverse scattering length passes through zero 
(via the tuning of the Feshbach resonance). We shall call 
it the halo dimer, in analogy with somewhat similar halo 
nuclei 0. A halo dimer is spatially extended and has 
only a small overlap with the closed-channel Feshbach 
resonance. But as the closed-channel resonance is tuned 
to lower energies, its energy eventually becomes negative 
and can be identified with the pair state. We present 
below a simple model that shows how this two-step tran- 
sition from the continuum pair states to the molecular 
bound state takes place, deriving the relevant size and 
energy scales. We then discuss the experimental signa- 
tures of the halo dimer state. Finally, we assess the the- 
oretical approaches that have been applied to treat the 
many-particle system in the halo dimer domain. There 
are a number calculations in the literature that ignore 
the special character of these state, and we believe that 
such models lead to a number of falsifiable predictions. 

We consider here a simple model of the two-body 
physics use a two-channel atomic Hamiltonian of the form 
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where r is the atom-atom separation. We shall use the 
units H = m = 1 and here u(r) describes the open two- 
atom channel and v(r) the closed (molecular) channel. It 
is implied that while Vn(r) — > and Vi 2 (r) = V 2 i(r) — > 
when r — » oo, the closed channel potential tends to a 
positive constant V 22 (r) — > Uq > 0. Since we are inter- 
ested in very low energies, only the s-wave is considered 
here and we show only the radial parts of the wave func- 
tions. If one were to apply this model to a system of 6 Li 
atoms, in the open channel at magnetic fields near the 
Feshbach resonance the two valence electrons of the two 
6 Li atoms would be anti-parallel to the magnetic field 
(triplet state), while in the closed channel the two elec- 
tron spins would be anti-parallel to each other (singlet 
state). As in all previous references [a lii LLii liii LL2| we 
shall assume that in the closed channel there is only one 
state close to the two-atom threshold. If the wave func- 
tion of this state is <f>o (r) with an energy Kq with respect 
to the two-atom threshold, one can easily show that the 
two-channel problem simplifies somewhat and the wave 
functions read 



u{r) = Uo(r) + ( r 
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v{r)=<f> (r) , 2 2 . 

Hu(r)u (r) = k 2 u (r), 
H 22 (r)(j) (r) = affair), 



-Vv. 



y o / ~r 2 zs~, (2) 



(3) 



(4) 
(5) 



For the sake of simplicity, we shall assume that in the 
open channel the atoms are free (the same assumption 
was made in Ref. 0), that (f>o(r) — sin(7rr/r )y / 2/Vo 
for < r < ro and zero otherwise and that Vi 2 (r) = 
Vzi(r) = gS(r — ri), with n = r /2. Physically, the pa- 
rameter ro is of the order of the van der Waals length 
fo ~ (Cem/H 2 ) 1 / 4 and we shall consider only such ener- 
gies for which kr <C 1. After some simple manipulations 
one arrives at the following form of the two-atom wave 
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function in the open channel (using uo(r) = sin(fcr)) 

g 2 exp(ifcr>) sin(fcr<) 



u{r) = sin(Ax) + 



Kq — g 2 — ig 2 kr\ 



(6) 



where r> = max(r, r\) and r< = min(r, r\) and terms of 
order 0(k 2 ) have been neglected. One can then easily 
show that 

-i = ^, (7) 
a 5 ri 

u(r) = exp(i(5) sin k(r — a), if r > r%, (8) 

where tan S — —ka. In particular, exactly at the res- 
onance u(r) = icos(fcr) for r > r\. Far away from 
the Feshbach resonance, the open channel wave func- 
tion would be approximately equal to u(f) — sin(fcr) 
instead (for r > ro). By the nature of the problem 
at hand we have Ko^o = 0(1) and gro = 0(1). Only 
if kq and g are comparable in magnitude one can at- 
tain the regime when a = ±oo. For the experiments 
reported, the relative energy of the Feshbach resonance 
and thus the magnitude of Kg is controlled by a magnetic 
field. With fine tuning one can make the numerator in 
Eq. J7J very small, and in this way attain the regime 
\a\ 3> ro, even though all the parameters in the equation 
above are of order 1/ro or ro respectively. As an order 
of magnitude estimate for these constants one can use 
0(kq) = 0(.g 2 ) = 2fiBBom/h 2 , where m is the atomic 
mass, fiB is Bohr magneton and Bq is the value of the 
magnetic field where 1/a = 0. For both 6 Li and 40 K one 
readily obtains that O(n r ) = O(gr ) = 1. 

Because of the coupling between the two channels, 
there is now a pole of the scattering amplitude at k = 
i/a, as may be seen from Eq. JfJJl. In different terms, 
by means of a magnetic field one controls the logarith- 
mic derivative of the open channel wave function near 
the origin, more exactly near r = ro. In a finite density 
medium it might naively appear that an infinite value 
of the scattering length would be meaningless. How- 
ever, that should be interpreted rather as the logarithmic 
derivative of the open channel wave function at r = ro, 
namely dhxu(r)/dr\ r — ro = —1/a, or in more physical 
terms, as the relative momentum with which the two 
atoms emerge after interacting at short distances. When 
the two-atom system is far from the Feshbach resonance 
the typical relative momentum with which the two atoms 
emerge after interacting at distances smaller than ro is 
h/ ro. The special situation, which is achieved by bringing 
the two atoms exactly at the Feshbach resonance, is to 
insure that they emerge from the interaction region with 
an essentially vanishing relative momentum. In a certain 
sense that amounts to an ultimate further cooling of the 
relative atomic motion to its minimum and that is what 
makes the physics of atoms under these conditions par- 
ticularly exciting. In a sense, two-atom collisions at short 
distances do not bring in any momentum. 



TABLE I: Character of the condensate as a function of the in- 
verse scattering length a -1 in various in intervals, the approx- 
imate boundaries of these intervals being shown in the second 
row. The total electron spin and spin projection (S, Sz) along 
the magnetic field for various pairs are shown in the last row. 
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Let us put this two-atom system in a spherical cavity of 
radius R ^> ro. In principle one would have to specify the 
boundary condition for u(r) at r = R, which would lead 
to energy quantization. The specific nature of this energy 
quantization (e.g. Neumann vs Dirichlet boundary con- 
ditions) is qualitatively unimportant. By choosing the 
radius of this spherical cavity one can simulate a Fermi 
gas of various densities, with kp oc 1/R, and in particular 
choose the regime with kp\a\ 3> 1, when according to the 
authors of Ref. the two atoms should be with prob- 
ability essentially one in the small size molecular state. 
One has simply to estimate the ratio of the probability to 
find the atoms in one or another channel, and one easily 
finds that 
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drsm z (kr + $) = O(R). 
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Consequently, the relative probability to find the two 
atoms in the molecular state is of the order 0(ro/i?) = 
O(kpro) <C 1 for a dilute Fermi gas, nrg <C 1. 

One can continue the argument in the region of neg- 
ative detuning, when a > 0, and easily convince oneself 
that as long as o > ro, the relative probability to find 
the two atoms in the closed channel is small. That was 
discussed in Ref . 1 131 and shown in another explicit cal- 
culation in Ref . 14] . Only when the scattering amplitude 
becomes of the order of the radius of the interaction (van 
der Waals length) the probability to find two atoms in the 
closed channel becomes comparable with the probability 
to find the same atoms in the open channel. The various 
regimes of the coupling strength are shown in Table I. 

Next we discuss experimental observables to distin- 
guish between the two kinds of bound states, i) One 
signature is magnetic. The experiments are typically 
carried out in magnetic fields large compared to the hy- 
pcrfinc splitting. The continuum states are thus polar- 
ized with respect to the electron spin, with the nuclear 
spin mainly responsible for the second component of a 
two-component Fermi gas. Thus the electron spin wave 
function of the halo dimer is largely S=l. The closed 
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channel, on the other hand, is typically well described 
with a spin-singlet electron wave function. In the exper- 
iment of Ref. |5|, the spin of the bound pair was mea- 
sured (see the insert of Fig. 2), and it was indeed found 
that its value was large down to values of the magnetic 
field well below the value at which the scattering length 
changes sign. This result is in perfect agreement with 
microscopic understanding of the Feshbach resonance in 
6 Li [J^. ii) Another signal of the character of the pair 
state in the condensate is in the particle loss rate of the 
system. When the condensate is prepared with the atoms 
in the lowest hyperfine states, the only inelastic processes 
are three-body collisions in which one of the pairs goes 
into a deeply bound molecule. Such processes are more 
probable when three particles are simultaneously in each 
other's range of interaction. The observations of Ref. 51 , 
see Fig. 4, confirmed by a separate experiment Ref. 3j, 
in fact show that the loss rate becomes large under the 
same conditions that the bound pair develops a singlet 
character. This observation also agrees with the theo- 
retical expectation of Ref. [rn |. Hi) We mention one 
more piece of evidence, this one more indirect. As sev- 
eral experimental groups have demonstrated Q when the 
scattering length is positive (thus on the BEC side of the 
Feshbach resonance) and when na 3 < 1 the sizes of the 
atomic clouds agree very well with theory based on halo 
dimers, which predicts that the interaction between them 
can be characterized by a scattering length of magnitude 
0.6a 0,^3- H the pairs were in closed-channel molecu- 
lar state, their scattering length would have the order of 
magnitude of their size C(^o) independent of the atom- 
atom scattering length a. 

We now turn to theoretical models for the conden- 
sation of trapped fermions into molecules. There have 
been a number of recent calculations [p,|p, Hoi Hl| , based 
on the model of Timmermans et ai, |l2j that assume a 
direct transition between the fermions and a molecular 
bound state. In Ref. [s] it was claimed that a dilute 
atomic Fermi gas near a Feshbach resonance undergoes a 
crossover into diatomic molecules of relatively small sizes 
on the BCS side of the resonance, when the atomic scat- 
tering length is still negative and the molecular energy 
level lies in the two-atom continuum. The authors of 
Ref. Q describe the results of the recent experiment of 
Regal et al. Q, by the way of this conversion of atoms 
into diatomic molecules. Another recent preprint makes 
a similar claim, that exactly this process occurs on the 
BCS side of the Feshbach resonance Q, although this 
particular work deals with different properties of such 
systems. 

The statement we take issue with here was formulated 
extremely succinctly by the authors of Ref. pj. One 
considers a uniform dilute Fermi gas of number density 
n = k F /3ir 2 , near a Feshbach resonance at positive de- 
tuning (a < 0), when the energy of the molecular state is 
e m ~ h 2 /ma 2 > 0. With respect to the two-atom contin- 



uum this molecular state is unbound. Then if kp\a\ > 1 
a fraction of these atoms converts into molecules with a 
number density, see Ref. and also Refs. |jj El where 
a similar result is quoted, 
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Naively this results seems obvious, since as soon as the 
Fermi energy exceeds the (positive) molecular state en- 
ergy, the system can only lower its total energy by con- 
verting into molecules. These molecules, being bosons in 
character, can all have zero center of mass kinetic energy, 
and thus the Fermi energy is lowered until it becomes 
equal to e m ~ h 2 /ma 2 and further conversion of atoms 
into molecules becomes energetically impossible. Apply- 
ing this to the two-body problem, the transition would 
take place discontinuously as the inverse scattering length 
approaches zero from the BCS side. But as we saw pre- 
viously, the formation of halo dimers is a smooth process 
with the two-particle wave function changing smoothly as 
the inverse scattering length goes through zero. It is also 
abundantly clear that such models 0, M El El would 
predict the wrong magnitude for the magnetic moments 
of such Fermi clouds near the Feshbach resonance. If in- 
deed the BEC of small size molecules would occur in this 
regime, the prediction would be that the magnetic mo- 
ment of such pairs would be vanishing, while experiment 
clearly shows that not to be the case p|. 

There is still one possibility, when the mechanism 
suggested in Refs. 0, |j| El El could prove indeed 
correct and near a Feshbach resonance such a system 
would condense into molecules of relatively small size 
O(ro). This situation would occur if the coupling con- 
stant g becomes unnaturally small. The relative proba- 
bility that two atoms spend most of the time in a small 
size configuration, according to Eqs. ({§1 - HOfl would be 
oc l/(g 2 riR) kpr /(gro) 2 , which could become acci- 
dentally large, even though the system is still dilute, in 
the sense that fcp-ro -C 1. This would happen if the Fes- 
hbach resonance would be accidentally extremely close 
to the atom-atom threshold. Whether this mechanism 
can be indeed realized in experiments with dilute atomic 
Fermi clouds, is a question that remains to be investi- 
gated. 

Since the halo dimers occupy the same sector of the 
two-particle Hilbert space as the plane waves of the Fermi 
gas, it is reasonable to ask whether a theory can be found 
that does not introduce a closed channel resonance ex- 
plicitly. In fact, the BCS theory has the degrees of free- 
dom to describe the many-particle system over the full 
range of couplings we consider here [lRj . One may show 
[lflj l that the BCS theory contains the correct description 
of the gas of halo dimers in the limit that a <C 1/kp 
and positive. Namely, the BCS chemical potential /i is 
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related to the dimer energy Ed by 

/i = T = -w- (12) 

On the other hand, the BCS theory is quantitatively in- 
correct on at least two properties. The pairing gap A 
is reduced by polarization effects from the BCS value by 
about a factor of « 2 both in the weak 20] and strong 
coupling limit |2lj . Also, the dimer-dimer scattering am- 
plitude is predicted to be 2a in the BCS theory . while 
a more detailed treatment gives 0.6a flEl flT^ . Still, it 
is still possible to use the a BCS framework for an ef- 
fective theory, renormalizing the couplings to reproduce 
these interactions, in the framework of the Superfluid 
LDA (SLDA) H3. 
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